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Summary 

A computer code for a single-point thermodynamic anal- 
ysis of hydrogen/oxygen expander-cycle rocket engines is 
described in this report. Full, split, and dual expander cycles 
are included. Heat regeneration between the turbine exhaust 
and the pump exhaust is allowed. The combustion process is 
modeled as one of chemical equilibrium in either an infinite- 
area or a finite-area combustor. Gas composition in the nozzle 
during expansion may be either equilibrium or frozen. 

This report, which serves as a users guide for the computer 
code, describes the system, the analysis methodology, and the 
program input and output and contains sample calculations to 
show the effects of key variables such as nozzle area ratio and 
oxidizer-to-fuel mass ratio. For this analysis, ideal specific 
impulse is a function of only nozzle area ratio and oxidizer- 
to-fuel mass ratio and is independent of expander-cycle type, 
chamber pressure, regenerator effectiveness, and turbomachine 
efficiencies. These latter parameters primarily affect the pump 
exit pressure required to achieve the specified chamber pressure. 

A single-point analysis is useful for determining system 
performance from known component performance and for 
performing sensitivity analyses over narrow ranges of the cycle 
variables. In analyzing a given engine over a wide range of 
cycle variables, care must be taken to use appropriate varia- 
tions in component performance parameters. 

Introduction 

Expander-cycle engines are the primary propulsion can- 
didates for orbital transfer missions. The advantages of an 
expander-cycle engine are good specific impulse, engine 
simplicity, and relatively low weight (ref. 1). For a topping 
cycle, such as the expander cycle, however, the fuel-pump 
discharge pressure is much higher than that for a bleed cycle, 
such as the gas generator cycle, and this tends to limit the 
achievable combustion chamber pressure. 

In support of space propulsion studies being conducted 
at the NASA Lewis Research Center, a computer code for 
the single-point thermodynamic analysis of hydrogen/oxygen 
expander-cycle rocket engines was developed. This code is 
applicable to full, split, and dual expander cycles, which can 
include a regenerator. Vacuum specific impulse, cycle tem- 
peratures, and cycle pressures are computed. Real fluid 
properties and methodologies are used for all computations 


except combustion and expansion. Combustion and nozzle 
expansion are modeled with one-dimensional, ideal-gas 
calculations. The combustion process is one of chemical 
equilibrium in either an infinite-area or a finite-area chamber. 
Expansion in the nozzle may be either an equilibrium or a 
frozen chemical process. 

This report presents the cycle analysis methodology and also 
serves as a users guide for the computer code. The expander 
cycles are described and the cycle analysis equations are 
presented. Program input and output are described and sample 
cases are included. The sample calculations are intended to 
show the effects of some of the key variables on the cycle 
performance. 

System Description 

The computer code presented herein can be used to analyze 
the single-point performance of full, split, and dual expander 
rocket engine cycles. Each of these cycles can incorporate a 
regenerator on the fuel side to transfer heat from the turbine 
exit stream to the pump exit stream. The three cycles are 
described in this section. The format for the system sketches 
is that used by Pratt & Whitney in reference 2. 

The full expander cycle with a regenerator is shown in 
figure 1 . Liquid hydrogen is pumped to high pressure and then 
flows through the regenerator, where heat is transferred to 
it from the turbine exit stream. The hydrogen is further heated 
as it regeneratively cools the nozzle and the combustion 
chamber. Most of the heated hydrogen then expands in series 
through the oxygen-pump and hydrogen-pump drive turbines. 
Alternatively, the order of the flow through the turbines can 
be reversed or the flow can be in parallel. Some of the heated 
hydrogen from the combustion chamber cooling jacket is 
bypassed around the turbines (and the regenerator) to provide 
a means for controlling turbine power. After transferring heat 
to the pump exit stream in the regenerator, the turbine exit 
flow mixes with the bypass flow, passes through the fuel 
shutoff valve, and enters the combustion chamber. On the 
oxidizer side the liquid oxygen is pumped to the pressure 
required to match the hydrogen pressure at the combustor. 
It then flows through the oxidizer control valve, which controls 
the oxidizer-to-fuel (o/f) ratio, and enters the combustion 
chamber, where combustion occurs. The combustion gases 
then expand through the nozzle to high velocity to provide 
the thrust. 
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A split expander cycle is shown in figure 2. This cycle 
derives its name from the fact that the hydrogen flow is split 
in the fuel-pump system. Only part of the hydrogen is pumped 
to maximum pressure for use as the chamber/nozzle regener- 
ative coolant and the turbine drive fluid. This figure shows 


Liquid 

hydrogen 



parallel flow through the turbines, which is one of the available 
modeling options. The remaining hydrogen flows through the 
jacket bypass valve, which controls the flow split, and mixes 
with the turbine exhaust and bypass streams before entering 
the combustion chamber. A regenerator can be incorporated 
into this cycle in a similar manner as for the full expander 
cycle (fig. 1). 

A dual expander cycle is shown in figure 3. In this cycle 
each of the fluids is heated by regenerative cooling of the 
nozzle or the combustion chamber, or both, and used in a 
separate turbine to drive its own pump. Thus, the fluids can 
be kept separated before combustion without relying on the 
turbopump seals. A bypass valve controls the flow to each 
turbine. If desired, a regenerator can be incorporated on the 
fuel side of the cycle. 

The method of analysis is documented later in this report. 
Symbols used for the analysis are defined in appendix A, and 
the analysis method is presented in appendix B. 

Program Input and Output 

Input 

The input data consist of namelist INPUT, which is required 
for the cycle calculations, and THERMO data, REACTANTS 
data, and namelists INPT2 and RKTTNP, which are used for 
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the combustion and expansion calculations performed by the 
Gordon and McBride CEC code (refs. 3 and 4). The CEC 
input data are described in detail in references 3 and 4 and 
are described here only to the extent necessary to run this 
program without consulting the CEC references. Input 
for a sample case is presented in figure 4. The output corre- 
sponding to this sample input is presented and described in 
the next section. 

Namelist INPUT is read from unit 2. The INPUT variables 
along with descriptions, units, and any special instructions are 
presented in the following list: 

PI fuel-pump inlet pressure (saturation 

temperature is assumed), psi 
P12 oxidizer-pump inlet pressure (saturation 

temperature is assumed), psi 
PC combustion chamber injection pressure, psi 

TBP fuel-turbine bypass fraction 


OBP oxidizer-turbine bypass fraction (dual 

expander only) 

JBP cooling-jacket bypass fraction (split expander 

only) 

TC thrust correction factor, ratio of actual thrust 

to ideal thrust 

OF oxidizer-to-fuel flow ratio 

ETAFP fuel-pump efficiency 

ETAFT fuel-pump drive turbine efficiency 

ETALP oxidizer-pump efficiency 

ETALT oxidizer-pump drive turbine efficiency 

ETABP low-pressure, fuel-pump efficiency (split 

expander only) 

EPS regenerator effectiveness (input only for 

IREG - 1) 


THERMO 

298.15 1000 
H 

0.25000000E 01 
0.25474390E 05- 
0.00000000 

H02 

0 . 40173060E 01 
-0.11412445E 04 

- 0 . 95674952E-08 
H2 

0 . 30558123E 01 
-0 . 861 68476E 03- 
0.74997496E-08- 
H20 

0 . 26340654E 01 
-0 . 29876258E 05 
-0.48670871E-08 
H202 

0 . 4792B858E 01 
-0 . 18132195E 05- 

- 0 . 46279780E-08 
0 

0.25342961E 01- 
0 . 29231108E 05 
-0.32604921E-08 
OH 

0.28897814E 01 
0 . 38857042E 04 
-0 . 52133639E- 09 
02 

0 . 361221 39E 01 
-0 . 1 19781 51E 04 
-0.98189101E-08 
03 

0 . 54665239E 01 
0 . 15214096E 05 
-0.20634230E-08 
H20CL) 

0.286308D0E 02 
-0 . 38579539E 05 
-0.13653020E-05 
END 

REACTANTS 
H 2. 

0 2. 


.000 5000.000 

J 3/7711 1. 

0 . 00000000 
0 . 45989041E 00 
o . oooooooo 

J 9/78H 1.0 

0.22175883E-02- 
0.37846051E 01 
0 . 36597628E-11 * 
J 3/77H 2. 

0.59740400E-03- 
0 . 17Z07073E 01 
0 . 25203379E- 1 1 - 
J 3/79H 2.0 

0 . 31121899E- 02' 
0 . 70823873E 01 
0.15284144E-11- 
L 3/85H 2.0 

0 . 3630 0865E- 02- 
0 . 51306415E 00 
0 . 21458200E-11 ■ 
J 3/770 1. 

0.12478170E-04 
0 . 49628591E 01 
0 . 1 0152035E-1 1 
J 6/770 l.H 
0.10005879E-02 
0 . 55566427E 01 
0 . 41826974E- 13 
J 3/770 2. 

0 . 748531 66E-03 
0 . 36703307E 01 
0 . 33031 825E- 1 1 
J 6/610 3.0 

0 . 1 7326031 E- 02 
■ 0 . 3471 2616E 01 
0 . 2001 5595E- 11 
J 3/79H 2.0 

• 0 . 2026 0986E 00 
-0.11895Q46E 03 
0 . 9 1 326966E- 09 


0. 0. 0.G 

0 . 00000000 
0.25000000E 01 
0 . 25474390E 05- 
2. 0 . O.G 

0 . 57710171E-06 
0 . 359641 02E 01 
0 . 89333502E 03 

0. 0. 0.G 

0 . 16747471E- 08- 
0 . 29432327E 01 
0 . 976954 1 3E 03- 

1. 0. 0 . G 

■ 0 . 90278449E- 06 
0 . 41675564E 01- 
0 . 30289546E 05- 

2. 0. 0.G 

■0.111 36435E- 05 

0 . 34546633E 01 
0 . 17672328E 05 

0. 0. 0.G 

■ 0 . 12562724E- 07 
0 . 30309401 E 01 
0 . 29136526E 05 

1. 0. O.G 

■0 . 22048807E- 06 

0 . 38737300E 01 
0 . 35802348E 04 
0. 0. O.G 

-0.19820647E-06 
0 . 37837 1 35E 01 
'0.1 06381 07E 04 

0. 0 0.0 O.G 

- 0 . 72204889E- 06 

0.24660617E 01 
0 . 16059556E 05 

1. 0. 0.L 

0.78529479E-03 
0 . 2863080 0E 02 

-0 . 38579539E 05 


300.000 5000 

0 . 00000000 
0 . 00000000 
0 . 45989841 E 00 
300.000 5000 

0 .71372882E-10 


000 


1 . 00790R1 


0.00000000 
0.00000000 
0 .00000000 
.000 33.00670 

0 . 3645859 1 E-14 


Formats 
3A4 
3F \ 0.3 

3A4,6X,2A3,4(A2 ( F3 0), A1 ,2F 1 0 3 ,F I 3 5.12 
5(E 1 5 8). 15 


0 .52500748E-03 0 . 75118344E-05 
0 . 66372671E 01 0.00000000 

300.000 5000.000 2.01580 

0 . 21 247544E- 1 0 0 . 25195487E- 1 4 
0 .3481 55 09E- 02-0 . 777 1 381 9E- 05 
0 . 18186137E 01 0.00000000 

300.000 5000.000 18.01520 

0 . 12673054E-09-0 . 69164732E-14 
0 . 181 06868E- 02 0 . 59450873E-05 
0 . 73087997E 00 0.00000000 

300.000 5000.000 34.01460 

0 . 14868513E-09-Q . 6895851 IE- 14 
0 .55575930E-02 0 . 921 03738E- 06 
0 . 68402452E 01 -0 . 16394994E 05 
300.000 5000.000 15.99940R1 

0.69029862E-11-0 . 63797095E- 15 2 

0.22525853E-02 0 . 39824540E- 05 3 

0 . 26099342E 01 0.00000000 4 

300.000 5000.000 17.00730R1 

0.20191288E-10-0.39409831E-15 2 

0 .13393772E-02 0 . 16348351E-05 3 

0 . 34202406E 00 0.00000000 4 

300.000 5000.000 31.99880 1 

0 .33749008E-10-0 . 23907374E-14 2 

■ 0 . 30233634E- 02 0 . 99492751E- 05 3 

0.36416345E 01 0.00000000 4 

300.000 5000.000 47.99820R1 

0 . 1 37 21 66 0E- 09-0 . 96233828E-14 2 

0 . 91703209E-02-0 . 49698480E-05 3 

0 . 12172130E 02 0.00000000 4 

273.150 500.000 18.01520 1 

-0.13653020E-05 0 . 91 326966E- 09 2 

-0.20260986E 00 0 . 78529479E-03 3 

-0.11895046E 03 0.00000000 4 


100.00 - 100 . 00 L 270.35 F 

100.00 -2000.00 L 129.86 0 


NAMELISTS 
& I NplJT 

pi = 18 0,P12 = 16.0,PC = 1500 . ,TBP = 0.05,TC = 1 .00 , DPCREG=0 . 02 ,DPHREG=0 . 04 , 
ETAFP=0.65,ETAFT=0.85,ETALP=0.70,ETALT=0.85,EPS=0.4, 

DPFI=0 . 15,DPLI=0 . 35,DTRISE = 450 . 0 ,DPIOSS=0 .15, 

0F=6 . 5 , CASE=1 , IH0= 1 , IREG=2 , IDBG= 0 , 

SEND 

8INPT2 P=1500 . , MIX=6 . 5 SEND 

SRKTINP EQL* .TRUE . ,FR0Z=.TRUE. ,FAC=.TRUE. , ACAT=2 . 5 , SUPAR = 1000 . > SEND 
STOP 


3A4 

A2,F7 5,35X,F7 5.19X.A1 


Figure 4.— Sample input. 
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DPCREG 

DPHREG 

DPFT 

DPLI 

DPBI 

DPLOSS 

DTRISE 

DPCHM 

DPNOZ 

DTCHM 

DTNOZ 

CASE 

IREG 

IHO 


regenerator cold-side fractional pressure loss 
(input only for IREG = 1) 

regenerator hot-side fractional pressure loss 
(input only for IREG = I) 

fuel-stream (line/valve/injector) fractional 
pressure loss 

oxidizer-stream (linc/valve/injector) fractional 
pressure loss 

bypass-stream (line/valve/injector) fractional 
pressure loss (split expander only) 

nozzle/chamber cooling jacket fractional 
pressure loss (full and split expanders only) 

nozzle/chamber coolant temperature rise, D R 
(full and split expanders only) 

chamber cooling jacket (fuel) fractional 
pressure loss (dual expander only) 

nozzle cooling jacket (oxidizer) fractional 
pressure loss (dual expander only) 

chamber coolant (fuel) temperature rise, D R 
(dual expander only) 

nozzle coolant (oxidizer) temperature rise, °R 
(dual expander only) 

cycle-type indicator 

1— full expander cycle 

2 — split expander cycle 

3— dual expander cycle 
regenerator indicator 

1— no regenerator 

2— fuel-side regenerator 

turbine arrangement indicator (full and split 
expanders only) 

1— series with fuel turbine followed by 
oxidizer turbine 

2— series with oxidizer turbine followed by 
fuel turbine 


3— parallel 

IDBG debug output indicator 

0— no additional output 

1— values of pressure /?, absolute temperature 
t, and enthalpy h for each station printed 
during calculations for both unconverged and 
converged passes 

The CEC input data are read from unit 5 and consist of 
THERMO data, REACTANTS data, and namelists INPT2 and 
RKTTNP. The THERMO data are the least-squares, curve-fit 
coefficients of the thermodynamic data for the species H, 
H0 2 , H 2 , H 2 0 (vapor and liquid), H 2 0 2 , O, OH, 0 2 , and 0 3 . 
These are the species considered for the equilibrium 
calculations for the H 2 /0 2 system. These THERMO data are 


used exactly as presented in figure 4. Formats are included 
in the table. Each time the THERMO data are read from 
unit 5, they are written in binary format to unit 4 and sub- 
sequently read from unit 4 for use in the program. Therefore, 
once the binary file exists, the THERMO data need not be 
included in the input data read from unit 5. 

The REACTANTS data identify each of the reactant species. 
These data are used exactly as presented in figure 4. Formats 
are included in the figure. 

The only variables to be included in namelist INPT2 are 

P must be set equal to value of PC in namelist 

INPUT 

MIX must be set equal to value of OF in namelist 
INPUT 

The variables that need be included in namelist RKTINP are 

EQL logical variable set to .TRUE, for equilibrium 
chemical composition in the rocket nozzle 

FROZ logical variable set to .TRUE, for frozen chemical 
composition in the rocket nozzle (defaults to 
.FALSE, if FAC - .TRUE.) 

FAC logical variable specifying combustion process 
.TRUE.— finite-area combustion, which includes 
total-pressure loss due to heat addition 
.FALSE.— infinite-area combustion 

ACAT ratio of combustor area to throat area (input only 
if FAC = .TRUE.) 

SUPAR rocket-nozzle supersonic area ratio (exit area to 
throat area) 

Output 

Program output consists of normal output, diagnostic output, 
and error messages. 

Normal output . — Normal output, which is written to unit 7, 
consists of a listing of input variable values, gas composition 
after both combustion and expansion, state conditions at all 
cycle points, and rocket performance. Figure 5 presents the 
normal output that corresponds to the sample input shown in 
figure 4. 

System configuration and input values are presented 
descriptively in the output (fig. 5). This particular case is for 
a full expander cycle with regenerator. The key parameters 
are a combustion chamber injection pressure of 1500 psia, an 
oxidizer-to-fuel ratio of 6.5, and a nozzle area ratio of 1000. 
All the component efficiencies and pressure losses are shown. 

This is specified to be a case with finite-area combustion 
and equilibrium expansion. Equilibrium gas compositions in 
terms of species mole fractions are shown for the combustor 
exit and the nozzle exit. The gas temperature at the nozzle 
exit is sufficiently low that the only species there are unburned 
hydrogen and steam. If the expansion were specified to be 
frozen, the nozzle exit composition would be the same as the 
combustion chamber composition. 
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ROCKET ENGINE CYCLE CODE 
FULL EXPANDER CYCLE WITH REGENERATOR 

FUEL PRESSURE = 18.0 PSIA 

OXYGEN PRESSURE - 16 ; ° £ S J A DCTrt 

INJECTION PRESSURE = 1500.0 PSIA 
THRUST CORRECTION FACTOR - 1.00 
FUEL TURBINE BYPASS RATIO = 0.05 
Q/F = 6.5 

sfW" o2 o T y?s iNE K Ys.r p 'ys' ator 

NOZZLE JACKET DELTA T = 450.0 R. 

S?lllllRi™5 E I»ELfA T ?." / fHiT , A»D CDLD SIDES) • 0 DA «ND O.OE 
FUEL INJECTOR (DELTA P)/P - 0.15 
LOX INJECTOR (DELTA P)/P = 0-35 


OPTIONS: 
CASE = 1 


IREG = 2 


I HO = 1 


EQUILIBRIUM MOLE FRACTIONS IN CHAMBER AND AT EXIT: 


H 
H02 
H2 
H20 
H202 
0 

OH 

02 


0.03274 
0 . 00010 
0.20194 
0.69675 
0 .00002 
0.00522 
0.05619 
0 . 00704 


0.00000 
0 .00000 
0.18102 
0.81898 
0 . 00000 
0 . 00000 
0 . 00000 
0.00000 


PRESSURE 
PSIA 

********** FUEL ********** 

FUEL PUMP INLET : 

FUEL HP PUMP EXIT : 

REGENERATOR COLD EXIT : 

COOLING JACKET EXIT : 

FUEL TURBINE EXIT : 

LOX TURBINE EXIT : 

REGENERATOR HOT EXIT : 

TURBINE BYPASS : 

MIXER EXIT : 

FUEL INJECTOR EXIT : 

******** OXIDIZER ******** 

LOX PUMP INLET : 

LOX PUMP EXIT : ?307.69 

LOX INJECTOR EXIT : 1500.00 

***** CHAMBER/NOZZLE ***** 

CHAMBER STAGNATION, EQL : 1452.30 

NOZZLE EXIT STATIC, EQL : 0.05 

EQUILIBRIUM PERFORMANCE PARAMETERS: 

NOZZLE AREA RATIO = 1000. 

CHARACTERISTIC VELOCITY = 7476. FT/SEC. 

NOZZLE EXIT MACH HUMBER = 6.912 . 

CHAMBER MACH NUMBER = 0.246 . 

VACUUM SPECIFIC IMPULSE = 493.4 LBF-SEC/LBM 


IB. 00 
3322.14 
3255.70 
2767.34 
2018.09 
1838.29 
1764.76 
2767.34 
1764.76 
1500.05 


CYCLE 

VALUES 

HORSEPOWER 

'TEMPERATURE 

FLOW 

DEGREES R 

LB/SEC 


37.75 

1 .000 

278.1 

79.62 

1.000 

307.58 

1.000 


757.58 

1.000 

278.1 

703.24 

0.950 

687.73 

0.950 

77.3 

436.78 

0.950 


757.53 

0.050 


452.33 

1 . 0G0 


453.17 

1 .000 


163.80 

6.500 

77.8 

175.54 

6.500 

178.72 

6479.23 

1191.71 

6.500 



Figure 5.— Sample output. 


Temperatures, pressures, and flow rates are shown for all 
cycle calculation stations. Flow rates are based on 1 lb/sec 
of fuel, since size effects are not considered in this single-point 
thermodynamic analysis. Also shown are turbomachinery 
power requirements. For this case, a fuel-pump exit pressure 


of 3322 psia is required to provide the specified combustion 
chamber injection pressure of 1500 psia. The finite-area 
combustion, which occurred at a chamber Mach number of 
0.246, resulted in a reduction in chamber stagnation pressure 
from 1500 to 1452 psia. Rocket performance is presented in 
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Pll AND P COMBUSTION ARE 
Pll AND P COMBUSTION ARE 
Pll AND P COMBUSTION ARE 
P1I AND P COMBUSTION ARE 
Pll AND P COMBUSTION ARE 


1632.03418 

1577.80420 

1521.67896 

1463.55029 

1500.04761 


1500 .00000 

1500.00000 

1500 . 00000 

1500 . 00000 

1500.00000 


Figure 6.— Chamber pressure convergence. 


Pll AND P COMBUSTION ARE : 1632.03418 

^? L -P. U M p _PRESSURE RATIO IS 208.333328' 


1500.00000 


FUEL PUMP EXIT . 
NOZZLE JACKET ENTR 
NOZZLE JACKET EXIT 
FUEL TURBINE INLET . 
POW1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
P0W1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
FUEL TURBINE EXIT : 
LOX TURBINE INLET : 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
LOX TURBINE EXIT : 
NOZZLE JACKET ENTR 
NOZZLE JACKET EXIT 
FUEL TURBINE INLET . 
POW1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
FUEL TURBINE EXIT ; 
LOX TURBINE INLET : 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
LOX TURBINE EXIT : 
NOZZLE JACKET ENTR 
NOZZLE JACKET EXIT 
FUEL TURBINE INLET . 
P0W1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
P0W1 AND DHTURB1 ARE 
P0W1 AND DHTURB1 ARE 
FUEL TURBINE EXIT : 
LOX TURBINE INLET : 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
LOX TURBINE EXIT : 
NOZZLE JACKET ENTR 
NOZZLE JACKET EXIT 
FUEL TURBINE INLET . 
POW1 AND DHTURB1 ARE 
POW1 AND DHTURB1 ARE 
P0W1 AND DHTURB1 ARE 
P0W1 AND DHTURB1 ARE 
FUEL TURBINE EXIT : 
LOX TURBINE INLET : 
P0W2 AND DHTURB2 ARE 
POW2 AND DHTURB2 ARE 
P0W2 AND DHTURB2 ARE 
LOX TURBINE EXIT 


3569.99976 
3498.59961 
2973.80957 
2973.80957 
210.422180 
210.422180 
210.422180 
210.422180 
2122.30225 
2122.30225 
55.0481567 
55.0481567 
55.0481567 
1932.98535 
3498.59961 
2973.80957 
2973.80957 
210.422180 
210.422180 
210.422180 
210.422180 
2122.64673 
2122.64673 
55.0481567 
55.0481567 
55.0481567 
55.0481567 
1933.39917 
3498.59961 
2973.80957 
2973.80957 
210.422180 
210.422180 
210.422180 
210.422180 
2122.72705 
2122.72705 
55.0481567 
55.0481567 
55.0481567 
1933.47485 
3498.59961 
2973.80957 
2973.80957 
210.422180 
210.422180 
210.422180 
210.422180 
2122.79712 
2122.79712 
55.0481567 
55.0481567 
55.0481567 
1933.58374 


Pll AND P COMBUSTION ARE : 1577.80420 

FUEL PUMP PRESSURE RATIO IS 198.333328 


82.6108398 

307.661865 

757.661865 
757.661865 

223.490814 " 
209.631226 
210.469269 
210.423935 __ 
699.821533 
699.821533 
56.9240112 
54.9870758 
55.0491791 
684.376953 

308.018555 

758.018555 
758 . 018555 

223.496719 
209 . 625900 
210.471451 
210.418610 
700.178467 
700.178467 
56.9275665 
54.9809570 
55.0568542 
55.0477905 
684.733398 

308.133545 

758.133545 
758.133545 

223.490814 
209 . 626694 
210.470078 
210 . 424728 
700.291016 
700.291016 
56.9182892 
54.9898224 
55.0497742 
684.844971 

308.169434 

758.169434 
758.169434 

223.495529 
209.635361 
210.464752 
210.423935 
700.326416 
700.326416 
56.9277496 
54.9912109 
55.0468140 
684.881836 


103.625839 _ 
953.772949 
2624.23120 
2624.23120 


2402.73218 

2402.73218 


2344.78564 

955.266357 

2625.48853 


2625.48853 


2403.99512 

2403.99512 


2346.05005 

955.746826 

2625.89307 

2625.89307 


2404 . 39331 
2404.39331 


2346 . 44604 
955.897217 
2626.02002 
2626.02002 


2404.52100 

2404.52100 


1500.00000 


2346.57690 


Turbine iterations 
to convergence 


Regenerator 
iteration loop 


Fuel-pump exit 
pressure 
iteration loop 


Figure 7.— Detailed diagnostic output. 


terms of characteristic velocity and vacuum specific impulse, 
which are 7476 ft/sec and 493 sec, respectively, for this case. 
Also shown is the nozzle exit Mach number, which is 6,9. 

Diagnostic output,— Diagnostic output is written to unit 6 
in two forms. A tracking of the convergence of combustion 
chamber injection pressure to the required input value is always 
printed as shown in figure 6. If the input variable IDBG = 1 , 
a longer diagnostic output, as shown In figure 7, tracking the 
entire computation is printed. Pressures, temperatures, and 


enthalpies are shown for the various calculation stations, and 
required-versus-calculated powers are shown for the turbines. 
Figure 7 shows one fuel-pump exit pressure iteration loop. 
The regenerator iteration loop is indicated as are the turbine 
iterations to convergence. 

Error messages .— The program contains three output 
messages that indicate abnormal operation due to either the 
nonexistence of a solution satisfying the specified input require- 
ments or the inability of a convergence algorithm to find a 
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solution. These messages are presented in this section, and 
their causes are discussed. In general, when one of the 
following messages appears, the program input should be 
checked for errors: 

(1) TURBINE N NOT CONVERGED!— The value of N 
is either l or 2 depending on which turbine in the sequence 
is not converged. This message is caused either by insufficient 
energy in the turbine drive fluid or by a failure in the conver- 
gence scheme. Corrective action includes reducing input chamber 
pressure or increasing turbine drive fluid energy, or both. 

(2) NO SOLUTION FOR CHAMBER PRESSURE, MAX 
PC = XXXXX.— This message is caused by the maximum 
achievable chamber pressure (shown in error message) being 
less than the chamber pressure specified by program input. 
If the input design parameters are correct, the input chamber 
pressure must be reduced. 

(3) CHAMBER PRESSURE SOLUTION NOT FOUND 
IN 100 ITERS— This message is caused by the program 
making 100 iterations without either finding a pump exit 
pressure that yields the required chamber pressure or 
determining that none exists. This situation has not occurred 
and there is no obvious reason for it to occur, but the possibility 
does exist. 


Sample Calculations 

Sample calculations were made to illustrate the effects 
of some key system parameters on cycle performance and 
operating conditions. The effects of nozzle area ratio and 
oxidizer-to-fuel mass ratio on ideal vacuum specific impulse 
were explored along with the effects of chamber pressure and 
regenerator effectiveness on pump exit pressure. 

The effect of nozzle area ratio (i.e., the ratio of exit area 
to throat area) on ideal vacuum specific impulse is presented 
in figure 8. Ideal specific impulse increased very rapidly with 
increasing area ratio up to an area ratio of about 200, after 
which the increase tapered off. Going from an area ratio of 
1000 to one of 2000 resulted in a 1 -percent increase in specific 
impulse. 

The effect of oxidizer-to-fuel mass ratio on ideal vacuum 
specific impulse is presented in figure 9. Maximum specific 
impulse occurred at an o/f ratio of about 6.5. Although 
combustion temperature was maximum at the stoichiometric 
o/f ratio of 8, the increase in reaction-product molecular 
weight with increasing o/f ratio caused the exit velocity, and 
therefore the specific impulse, to reach a maximum at an o/f 
ratio somewhat less than stoichiometric. 

Figures 8 and 9 were discussed without reference to a 
specific cycle (full, split, or dual) because those results were 
independent of cycle type. For this analysis, ideal specific 
impulse is a function of only nozzle area ratio and o/f ratio 
and is independent of cycle type, chamber pressure, regenerator 
effectiveness, and turbomachine efficiencies. These latter 



Figure 8.— Effect of nozzle area ratio on ideal vacuum specific impulse. 
Chamber pressure, 1500 psia; oxidizer-to-fuel ratio, 6.5. 



Oxidizer-to-fuel mass flow ratio 

Figure 9.— Effect of oxidizer-to-fuel mass flow ratio on ideal vacuum specific 
impulse. Chamber pressure, 1500 psia; nozzle area ralio, 1000. 

parameters primarily affect pump exit pressure, as illustrated 
in the next figure. 

The effects of chamber pressure and regenerator 
effectiveness on fuel-pump exit pressure for the full expander 
cycle are presented in figure 10. Increasing the chamber 
pressure obviously required higher pump exit pressures. 
Increasing the regenerator effectiveness allowed a lower pump 
exit pressure for a specified chamber pressure or, conversely, 
a higher chamber pressure for a given pump exit pressure. 
For example, a chamber pressure of 1500 psia required a 
3860-psia pump exit pressure with no regeneration but only 
a 2800-psia pump exit pressure with regeneration of 0.6 
effectiveness. Turbine inlet temperature increased with 
regeneration, thus reducing the turbine pressure drop required 
to provide a given power output. 
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Regenerator 

effectiveness, 



Chamber pressure, psia 

Figure 10.— Effects of chamber pressure and regenerator effectiveness on fuel- 
pump exit pressure (full expander cycle). Oxidizer-to-fuci ratio, 6.5; nozzle 
area ratio, 1000. 

Concluding Remarks 

In a single-point thermodynamic analysis all component 
performance parameters (i.e., turbomachine efficiencies, heat 
transfer effectiveness, and fractional pressure losses) are input 
as single values. Such an analysis, with the component param- 
eters remaining constant, is used to obtain system performance 
at a single point where the component performances are 
known. It is also useful for performing sensitivity analyses 
over narrow ranges of the cycle variables where component 
performance would remain essentially constant. 

Over a wide range of cycle variables and system size, it 
would be expected that component performance would vary 


as the component design and off-design conditions (e.g., 
chamber pressure, pressure ratio) varied. Therefore, for 
accurate analysis of an engine over a wide range of conditions, 
variations in the component performance parameters, as 
determined externally to this code, should be included as part 
of the program input. Future enhancements of this code could 
include incorporation of algorithms to correlate component 
performance with the cycle variables. 

Summary of Results 

This report presents a single-point thermodynamic code for 
analyzing hydrogen/oxygen expander-cycle rocket engines. 
The code is capable of analyzing full, split, and dual expander 
cycles. Heat regeneration between the turbine exhaust and the 
pump exhaust is allowed. The combustion process is modeled 
as one of chemical equilibrium in either an infinite-area or 
a finite-area combustor. Expansion in the nozzle can be either 
equilibrium or frozen. 

This report, which serves as a users guide for the computer 
code, describes the system, the analysis methodology, and the 
input and output. Sample calculations are included to show 
the effects of some key variables. These sample calculations 
showed 

1 . That ideal specific impulse for this analysis is a function 
of only nozzle area ratio and oxidizer-to-fuel mass ratio and 
is independent of expander-cycle type, chamber pressure, 
regenerator effectiveness, and turbomachine efficiencies 

2. That cycle type, chamber pressure, regenerator effective- 
ness, and turbomachine efficiencies primarily affect the pump 
exit pressure required to achieve the specified chamber pressure 

This single-point analysis is useful for determining system 
performance from known component performance and for 
performing sensitivity analyses over narrow ranges of cycle 
conditions. For the analysis of an engine over a wide range 
of cycle variables, appropriate variations in component 
performance should be included. 
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H GASPLUS function for enthalpy, Btu/lb 

h enthalpy, Btu/lb 

off oxidizer-to-fuel mass flow ratio 
P power, Btu/sec 

p pressure, psi 

S GASPLUS function for entropy, Btu/(lb)(°R) 

s entropy, Btu/(lb)(°R) 

T GASPLUS function for temperature, °R 

t absolute temperature, °R 

w mass flow rate, lb/sec 

x mass fraction 

e regenerator effectiveness 

7j efficiency 

Subscripts: 

b bypass 

CEC referring to CEC code 

/ fuel 

fcj fuel cooling jacket 


Appendix A 
Symbols 

fp fuel pump 

ft fuel turbine 

fib fuel turbine bypass 

id ideal 

inj injection 

jbp jacket bypass 

loss loss 

Ip low pressure 

o oxidizer 

ocj oxidizer cooling jacket 

op oxidizer pump 

ot oxidizer turbine 

otb oxidizer-turbine bypass 

rc regenerator cold side 

rh regenerator hot side 

sat saturated 

t turbine 

0-19 station identification (see figs. 1-3) 



Appendix B 
Analysis Method 


This analysis makes use of the GASPLUS properties 
program (ref. 5) available at NASA Lewis. The thermo- 
dynamic property functions obtained from GASPLUS and 
incorporated in the computer program are hereinafter 
expressed in the form 

y = F(z x ,z 2 ) ( 1 ) 

where y is the value of the function F t and Z\ and Zi are the 
independent variables. These are real-fluid properties and are 
used for all computations except those for combustion and 
nozzle expansion. The previously referenced ideal-gas equili- 


brium analysis code (CEC code) of Gordon and McBride (refs. 3 
and 4) is used for the combustion and expansion analyses. 

The analysis method for the full expander cycle, as shown 
in figure 1, is presented in full. Then, the modifications 
required for the split and dual expander cycles are presented. 
The calculation procedure for these cycles, particularly for 
the fuel system, is iterative, involving two or three levels of 
iteration depending on whether there is regeneration. For the 
oxidizer system, only the dual expander cycle requires such 
iteration. A calculation logic diagram illustrating the nature 
of this iterative procedure for the fuel system is shown in 
figure 1 1 . 



Figure 1 1 .—Fuel -system calculation logic. 
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Full Expander 

Starting on the oxidizer side, oxygen enters the system at 
input pressure p \ 2 and saturation temperature 1 \ 2 , where 

*12 - ^s»t(P 1 2 ) (2) 

Oxidizer enthalpy, entropy, and flow rate are then given by 


t\ = Tjipi) 

(13) 

Fuel enthalpy and entropy are then given by 


/?, = H(t u p{) 

(14) 

and 


s\ = 5(r,,pi) 

05) 


h\ 2 = H{t n ,p n ) 

(3) 

J i2 = S(t\2<Pn) 

(4) 


w 12 * (o/f)w | (5) 

where o//is the input oxidizer-to-fuel mass flow ratio and w, 
is the fuel mass flow rate, whose value is set to unity for this 
analysis. 

In the oxidizer pump the oxygen is pumped to pressure /? l3 , 
where 


The hydrogen is pumped to pressure p 2 in the fuel pump. A 
converged value of p 2 must be determined iteratively on the 
basis of the calculated hydrogen pressure after injection 
matching the input value of combustion chamber injection 
pressure p inj . On the basis of an initially estimated value of 
p 2 , the fuel-pump exit conditions are determined as 


^2, id “ H(S\,p 2 ) 

(16) 

A j ^2, id “ h\ 

Ah jp = 

(17) 

njj, 


h 3 — + &hjj, 

(18) 


P 13 * 



( 6 ) 


h — T(h 2 ,p 2 ) (19) 

The thermal performance of the regenerator is specified by 
its effectiveness, which is an input and is defined as 


The combustion chamber injection pressure p m} and the 
oxygen-stream pressure loss ratio Ap 0 \vdP\i> which represents 
the sum of the line, valve, and injector pressure losses, are both 
input. Pump exit conditions are then determined from 


*!3.id = 

(7) 

, , *13, id — ^12 

Ah op = — 

(8) 

Vop 


h\j ~ *12 + Ah » P 

(9) 

fi3 = T(h ]3 ,p i3 ) 

(10) 


__ W 2^1 ~~ ^2) _ ^19 ^9 ^20) 

Wi9(A|9 — h 2 ) h\q — h 2 

where the hot-side flow rate w, 9 is less than the cold-side 
flow rate vv 2 as a result of the fuel-turbine bypass ratio xj^,, 
which is input. 

W|9 = (1 ~Xftb)w 2 < 21 ) 

Since the value of /i 3 depends on the value of h i9 , which is 
unknown at this point, these values must be determined 
iteratively. With an initial estimate for h ]9 , we get 


Assuming the pressure loss Ap <}>] oss to occur at constant 
enthalpy, the oxygen enthalpy and temperature after injection 
into the combustion chamber at pressure p m} are 


*3 


= h 2 + e 



( 22 ) 


ho,ir\} h\i 


(U) 


and 


^j.inj = T(h ti jnj^Pinj) (12) 

On the fuel side, hydrogen enters the system with flow rate 
w, at input pressure p, and saturation temperature where 


p y =p 2 [ 1 - 


A Prc 

Pi 


(23) 


where the regenerator cold-side pressure loss ratio A p n is 
input. Then 


*3 = 


(24) 


m 



Note that with no regenerator (i.e., e = 0 and Ap r( .!p 2 = 0), 
the variable values at station 3 are equal to those at station 2. 

The hydrogen is then used to regeneratively cool the nozzle 
and the combustion chamber. Fuel cooling-jacket temperature 
rise and pressure loss ratio A Pfijlpi are input, and the 
cooling-jacket exit conditions are given by 


U = h + A t fi .j 

( 25 ) 

-4 1 - Afc) 

(26) 

V Pi ) 


h 4 = H(t 4 ,p 4 ) 

(27) 

s 4 = S(t 4 ,p 4 ) 

(28) 


Some of the hydrogen leaving the cooling jacket is bypassed 
around the turbines, and the hydrogen flow rate entering the 
turbine system is the same as that previously determined for 
the regenerator hot side. That is, 

h’5 = W|9 (29) 

The state conditions at station 5 are the same as those at 
station 4. 

There are three options for specifying the turbine system 
configuration: (1) series with oxidizer-pump drive turbine 

first, as shown in figure 1, (2) series with fuel-pump drive 
turbine first, and (3) parallel, as shown in figure 2. For any 
of these configurations the turbine power requirements are 


Pfi = W’2 a hjp (30) 

P,„ = w l2 Ah op (31) 

and the turbine specific works are 

P* 

A h fi = (32) 

w ft 

A ft* = — (33) 

Wot 

For the series arrangements 

= W,„ = H’5 (34) 


For the parallel arrangement the flow w 5 is split between the 
two turbines on the basis of equal pressure ratios (i.e., equal 
ideal specific works). This yields 

A/iij^ = A h idj „ (35) 


_ r<>, 

WfiVfl W i>/Vnr 

(36) 

W J1 _ PjiVo, 
W Of Po'Vfi 

(37) 

and since 


w 5 - w fi + w ot 

(38) 

we get 


PfiVot 

Wjj — M’ 5 

P o,Vfi + 

(39) 

and 


= w 5 - w fi 

(40) 

With the required specific work now determined for each of 

the turbines, the solution procedure is the 

same for each of 

the two turbines in any of the three arrangements. An iteration 

to determine turbine pressure ratio is 

required, since 

GASPLUS does not provide a means for calculating pressure 
with enthalpy and entropy given as the independent variables. 

Using the turbine arrangement in figure 1 

as an example, a 

pressure ratio p 6 fp 5 across the oxidizer turbine is assumed and 


(41) 

h 6M = H(s 5 ,p 6 ) 

(42) 

A/f„, Votih 5 ^6, id) 

(43) 

The turbine pressure ratio is varied and equations (41) to (43) 
are repeated until the calculated value of turbine specific work 
A h ol from equation (43) matches the required value from 

equation (33). Then, 


<1 

I 

II 

■*£! 

(44) 

'6 = T(ht>'Pb) 

(45) 

s b = S(! b ,p b ) 

(46) 


The fuel-turbine inlet conditions at station 7 are the same as 
the oxidizer-turbine exit conditions at station 6. A similar 
procedure for the fuel turbine then yields /? 8 , /* H , and ; 8 . In 
the parallel arrangement the inlet conditions for the fuel turbine 
are those at station 5 rather than those at station 7. 

Where the turbines are in series, the regenerator hot-side 
inlet conditions (station 19) are equal to the turbine-system 
exit conditions (station 8). With the turbines in parallel the 
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exit flows are mixed such that 


and 


P 19 — Pb ~ P* 


(47) 


r /inj = T(/i /inj ,/?inj) ( 59 > 


h ]9 = W ° thb + Wj ~ (48) 

r l9 = T(hw,pw) ( 49 ) 

The calculated value of /i 19 is compared with the estimate 
used for the regenerator cold-side evaluation in equation (22), 
and equations (22) to (48) are iterated to convergence. 
Regenerator hot-side exit conditions can then be determined 
from the following: 


. A/> W| \ 

P? =PM 1 ) 

\ P 19/ 

(50) 

ht) — h\g e (/1 19 — h 2 ) 

(51) 

t 9 = T()h h p<>) 

(52) 

The regenerator hot-side exit flow is then 
turbine bypass flow 

mixed with the 

H'io = "’5 + Xj,h w 2 = m ’2 

(53) 

P\o = Pv 

(54) 

, W’,/!y + Xjj h W 2 h 4 

«I0 “ 

H’io 

(55) 

cx. 

0 

11 

0 

(56) 


The fuel-stream pressure loss ratio AP/.iosiZ/’io? which is an 
input representing the sum of the line, valve, and injector 
pressure losses, is assumed to occur between the mixer and 
the combustion chamber. A combustion chamber injection 
pressure for the fuel is then determined as 


The combustion and nozzle expansion analyses are per- 
formed by using the program CEC, described in references 3 
and 4, as a subroutine. All the capabilities of the CEC program 
are available to this analysis (i.e., infinite-area or finite-area 
combustion, equilibrium or frozen expansion). The analysis 
methods are presented in the references and are not described 
herein. The CEC program uses an ideal-gas enthalpy reference 
temperature of 536.67 °R, which differs from the GASPLUS 
reference temperatures. Therefore, the CEC reactant 
enthalpies must be corrected. Furthermore, the enthalpy gained 
by the hydrogen during regenerative cooling of the nozzle and 
the combustion chamber must be accounted for. This is 
modeled herein by subtracting that enthalpy gain from the 
hydrogen enthalpy. Therefore, the CEC reactant enthalpies are 

h ft cEc = h inj - 77(536.67 °R, 14.696 psia) - (/t 4 ~ *3) 

(60) 

ho f CEC = h oM - 7/(536.67 °R, 14.696 psia) (61) 

It is recognized that the CEC program thermodynamic data 
are for normal hydrogen whereas the rocket-engine fuel is 
parahydrogen. Over the temperature range of combustion and 
expansion, the difference in the computed value of ideal 
specific impulse would be very small, if not insignificant. The 
ideal specific impulse computed by the CEC program is then 
multiplied by the thrust correction factor to yield the actual 
specific impulse. 

Split Expander 

The split expander cycle, shown in figure 2, differs from 
the full expander cycle in that the hydrogen flow is split in 
the fuel-pump system with part of it bypassing directly to the 
combustion chamber. In the low-pressure fuel pump all of the 
hydrogen is pumped to intermediate pressure p 15, where 


PfM 



Ag/Uoss \ 

P10 / 


(57) 


P 15 = 


7* inj 



(62) 


The fuel-pump pressure p 2 is varied and equations (16) to 
(57) are repeated until the calculated value of from 
equation (57) matches the input value of p in This iteration 
is controlled by the CONTIN subroutine described in 
reference 6. Then 


The combustion chamber injection pressure p mj and the 
hydrogen-bypass-stream pressure loss ratio Ap bM Jp^, which 
represents the sum of the line, valve, and injector pressure 
losses, are both inputs. Pump exit conditions are 


fy, inj *10 


(58) 


* 15, id — H(S[,p 15 ) 


(63) 


13 


(73) 


A/l _ ^IS.id - h \ 

^fp'lp 

(64) 

A hfpjp 

(65) 

= T(h [5 ,p ]5 ) 

(66) 

s \5 = S(f|5,p| 5 ) 

(67) 

The flow leaving the low-pressure fuel pump is split as 
specified by the input bypass fraction jc jhp . An amount, 
termed jacket bypass flow u j 15 , where 

vv !5 “ *jbp w 1 

(68) 


P 16 ~ P 13 



^16 “ H(hb>P\t) 

(74) 

s \6 = ^ (t\6>P\b) 

(75) 


Some of the oxygen leaving the cooling jacket is bypassed 
around the turbine, and the flow through the turbine is 

w 17 “ w I2 (l -x ofh ) (76) 

where x olb is the input oxidizer-turbine bypass ratio. The 
oxidizer-turbine required specific work is 


bypasses the regenerator, the nozzle/chamber cooling jacket, 
and the turbines and rejoins the main flow downstream in the 
mixer. The flow rate and enthalpy leaving the mixer are 


A h m * 


Wj7 


W I2^ 

VV |7 


(77) 


W,o = W 5 + X fih W 2 + W|5 

(69) 

, w s h 9 + x fih w 2 h 4 + W| 5 /t|5 

"10 

(70) 


W-’IO 


For this cycle the corrected fuel enthalpy for the CEC program is 

hj,c ec = inj - 7/(536.67 °R, 14.696 psia) 

“0 ~Xjhp){h A -h 3 ) (71) 

The analysis of the high-pressure fuel pump is the same as 
that for the full expander cycle except that the pump inlet 
conditions are those for station 15 rather than station L 


and a turbine pressure ratio P\j/p\ b is estimated. Then, 


<7 

•^L> 

sO 

ft. 

11 

r-> 

ft. 

(78) 

^17, id ^ H( s \b>P\l) 

(79) 

A hot = Voi(fo 16 ^I7,id) 

(80) 


the turbine pressure ratio is varied, and equations (78) to (80) 
are repeated until the calculated turbine specific work A h of 
from equation (80) matches the required value from 
equation (77). The turbine-exit enthalpy and temperature are 
then given by 


Dual Expander 

The fuel-system analysis for the dual expander cycle, shown 
in figure 3, is the same as that for the full expander cycle except 
for the simplification of having only one turbine in the system. 
Therefore, the variable values at stations 6 to 8 and 19 are 
identical. 

The oxidizer system is similar to the fuel system except that 
there is no regenerator. The oxygen-pump exit pressure /?j 3 
must be determined iteratively in the same manner as the 
hydrogen-pump exit pressure. An estimated value of /?, 3 is 
used to determine oxidizer-pump exit conditions at station 13 
in the same manner as for the full expander cycle. The oxygen 
is then used to regeneratively cool the nozzle. The cooling- 
jacket temperature rise A t acJ and the pressure loss ratio 
&PacjtP\s are input, and the nozzle cooling-jacket exit 
conditions are given by 

r i6 = *13 + A t (K .j (72) 


^17 “■ ^16 “ A h of 

(81) 

hi m T(h n> p n ) 

(82) 

The turbine flow and the bypass flow then mix 


W| 8 = W I2 

(83) 

P\% Pm 

(84) 

L ^17*17 + X ofh W ]2 

"18 - 

^18 

(85) 

^18 = 7’( J ^l8tPl8) 

(86) 


The oxidizer-stream pressure loss ratio A/? fjJoss /p , 8 is assumed 
to occur between the mixer and the combustion chamber. A 
combustion chamber injection pressure for the oxygen is then 
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determined as 


^.inj — 


(89) 




(87) and the corrected oxygen enthalpy for use in the CEC program is 


and compared with the input value of combustion chamber 
injection pressure p iny The oxidizer-pump exit pressure p u 
is varied until the calculated and input injection pressures 
match. Then, 


K,cec = Kms ~ H (536.67 °R, 14.696 psia) - (A, 6 - M 

(90) 


^.inj - /? I8 


( 88 ) 
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